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Abstract. We provide a fully quantum description of a mechanical oscillator in the presence 
of thermal environmental noise by means of a quantum Langevin formulation based on 
quantum stochastic calculus. The system dynamics is determined by symmetry requirements 
and equipartition at equilibrium, while the environment is described by quantum Bose fields 
in a suitable non-Fock representation which allows for the introduction of temperature. A 
generic spectral density of the environment can be described by introducing its state through 
a suitable P-representation. Including interaction of the mechanical oscillator with a cavity 
mode via radiation pressure we obtain a description of a simple optomechanical system in 
which, besides the Langevin equations for the system, one has the exact input-output relations 
for the quantum noises. The whole theory is valid at arbitrarily low temperature. This allows 
the exact calculation of the stationary value of the mean energy of the mechanical oscillator, 
as well as both homodyne and heterodyne spectra. The present analysis allows in particular to 
study possible cooling scenarios and to obtain the exact connection between observed spectra 
and fluctuation spectra of the position of the mechanical oscillator. 
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1. Introduction 

Optomechanical systems in the quantum regime are very important for quantum information 
processing and for testing fundamental issues of quantum mechanics Their theoretical 

analysis therefore calls for a first principle description. In particular since the focus is on 
quantum effects, the theoretical models must be fully consistent with quantum mechanics. 
Actually the correct quantum description of a mesoscopic mechanical oscillator and of the 
thermal noise affecting it is not a trivial task, and there is not a unique accepted model for 
them GM- 

The first aim of this paper is therefore to obtain an accurate quantum mechanical 
description of a mechanical oscillator taken to be part of an opto-mechanical device. 
The oscillator cannot be considered as a Brownian particle, but rather as a mesoscopic 
mechanical system, say a movable mirror mounted on a vibrating structure. Dissipative 
effects are essentially due to the interaction with phonons. Our strategy will be to introduce 
reasonable physical requirements leading to a master equation in Lindblad form, valid for 
any temperature of the thermal bath. We then translate these results into quantum Langevin 
equations and we show how to obtain a suitable non-Markovian generalization at this level 
of description. Relying on these results we can consider the description of the simplest 
optomechanical system, that is a moving mirror interacting with an electromagnetic mode 





















CONTENTS 


3 


in a cavity via radiation pressure p]|5]-[7||25|. Again a suitable analysis of the composite 
system and of the monitoring of the emitted light calls for a consistent quantum description. 
We shall obtain this result by the use of quantum Langevin equations, directly deducing them 
from a unitary dynamics, and exploiting the theory of measurements in continuous time. 

The paper is organized as follows. In Section]^ we determine the reduced dynamics 
of the mechanical oscillator. Here, the basic assumption is the use of a Markovian master 
equation with a quadratic generator and having a unique equilibrium state. Its structure is 
further determined by suitable symmetry requirements and by physical constraints on the 
behaviour of the mean values of position and momentum. In Section]^ we introduce the 
quantum Langevin equations for the mechanical oscillator alone; the whole presentation is 
based on the notions of quantum noise | |2^[Z7) and of input-output fields p8|j30l , as well as 
on the use of quantum stochastic calculus [ |3T|[32) . The Bose helds entering in the unitary 
dynamics play the role of phonon fields. By modifying their state without changing the time 
evolution operator it is possible to introduce non Markovian effects, namely a non-flat noise 
spectrum. The differences with respect to usual approaches are of relevance especially at low 
temperatures, where the zero-point fluctuations play an essential role. 

A quantum optomechanical system is studied in Section by using the quantum 
Langevin approach, within a fully consistent formulation valid at any temperature. Firstly, the 
typical effect of laser cooling is discussed. Then, the continuous monitoring of the emitted 


light is introduced in Sections 4.3.1 (homodyne detection) and 4.3.2 (heterodyne detection). 


The treatment is well based in the theory of measurements in continuous time. Detection 
of the emitted light is usually assumed to give a direct measurement of the fluctuations of 
the position of the mechanical component. We show that this is true, but only for not too 
low temperatures; at very small temperatures, interference terms are important and the direct 
connection with such fluctuations is lost. This leads to new predictions on the optical spectra 
at very low temperature. We hnally summarize and discuss our results in Section]^ 

2. Damped mechanical oscillator: the master equation approach 


As a first step towards the construction of models of optomechanical systems valid in 
the quantum regime at low temperatures, we consider the reduced dynamics of an open 
mechanical oscillator. A fully consistent quantum description of a massive nanomechanical 
component, kept at the simplest possible level, will be our basic building block in order to 
consider more complex dynamics. We therefore formulate in the first instance a Markovian 
description for the mechanical oscillator, which we build up relying on general physical 
constraints and symmetry requirements. 

A standard approach often considered in the literature is to derive the master equation 
for the harmonic oscillator from effective environmental models of bosonic oscillators | [T3][T4l 
20 ■ 2^ . However, previous work imiizi las shown that, while using careful approximations 
a positive Markovian dynamics can be obtained in this framework, the hnal results are valid 
only from medium to high temperatures of the thermal bath. To ask for a Markovian dynamics 
based on symmetry arguments allows to get simpler and more universal models, but again 
serious problems appear. The requirement that the equilibrium state should be the canonical 
thermal state determined by the standard Hamiltonian of a harmonic oscillator is known to 
be incompatible with positivity and translational invariance | [TT][3^ E)- This incompatibility 
induced some authors to renounce to translational invariance | |15|34| , or to accept non-positive 
dynamical equations and to give more relevance to obtaining time evolutions very close to the 
classical ones 113 14 A non positive dynamics can be satisfactory when the system is 
near the classical regime, but this approach becomes questionable when quantum effects are 
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searched for | |3^[J7) . We shall show that it is possible to maintain positivity and translational 
invariance by weakening the requests on the equilibrium state. The key point will be a weak 
formulation of energy equipartition at equilibrium. Our result is therefore to single out from 
the many proposals appearing in the literature a unique consistent dynamics in the Markov 
approximation. 


2.7. Physical constraints and symmetry requirements 


We formulate now our assumptions, starting from the existence of a well dehned positive 
Markovian dynamics, describing damping and translationally invariant apart from the 
harmonic potential. A weak equipartition condition and the existence of a unique stationary 
state in Gibbs form, as we shall see, will essentially fix the structure of the reduced dynamics. 
Assumption 1 (Positive Markovian dynamics with quadratic generator). The evolution of the 
statistical operator of the oscillator is governed by a Markovian master equation preserving 
the positivity of the states. The generator of the dynamics is at most quadratic in the position 
and momentum operators of the mechanical oscillator. 


The first assumption is to consider a time-homogeneous and linear time evolution. Such a 
dynamics can be expressed in the form 


^p(t)=£[p(f)], 


( 1 ) 


with L a suitable generator or Liouville operator, at most quadratic in the position and 
momentum operators of the mechanical oscillator q and p, so as to have at most a quadratic 
potential term and a friction effect proportional to the momentum of the mechanical oscillator. 
In the case of linear systems it is known that positivity and complete positivity of the dynamics 
are actually equivalent ini^ therefore according to the generator C must have the 

standard Lindblad structure. The most general quadratic Liouville operator is obtained in 
terms of two Lindblad operators in) 


Rj = {ujq + Vjp) , 


u,v G 


( 2 ) 


and a generic selfadjoint quadratic Hamiltonian for the mechanical system 

= y 9^ + ^ {<?,p} + y + fqq + fpP, 
where all the constants are taken to be real, so that C takes the form 

C[p] = -^ [77^, p] + ^ (r,pR] - ^ p}) • (3) 

7=1 ^ ^ 

We ask now to have a damped oscillator, but not an overdamped one. 


Assumption 2 (Damping). The “kinetic energy” term is non negative and the mean values of 
position and momentum decay to zero with an oscillating behaviour. 


Apart from the trivial requirement of a positive kinetic energy term, we further look for a 
dynamics describing the oscillating decay of the mean values of q and p to zero. This condition 
complies with the Markovian and quadratic approximations, which are expected to be good 
only for small damping. Denoting by {X)t the mean value of a quantum operator with the 
state p{t) solution of the master equation we have for position and momentum 

= hp{p)t + (y - Im {u\v)'j {q)t + fp, 
yy = -hq{q)t “ (y + Im(wk)) (p)t “ fq, 
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with (wit;) the inner product in C^. The eigenvalues of the associated dynamical matrix are 
— Im (u|t;) ± — hphq, so that in order to have an underdamped dynamics we need 

Im {u\v) > 0 and Ko/4 < hphq. In particular hp and hg have the same sign and are non¬ 
vanishing. Positivity of the kinetic energy leads to hp > 0 and therefore hg > 0. Then, we can 
write hp = 1/m and hg = the above inequality on kq becomes Kq < 417^. Finally, 

the vanishing of the equilibrium means imply fg = 0, fp = 0. Introducing the positive 
coefficients 

7 i„ = 2Im(M|u), Dqp=Re{u\v), Dgg = WvW"^ , Dpp = WnW"^ , 
the generator can be written in the form 

1 J~) T) 

mp] = - 2 ^ b, {p, p}] - ^ [9, p}] - b, b, p]] - ^ [p, [p, p]] 

- Dqp [p, [q, p]] - ^ ^ [q, {p, p}] - ^ ^ [p, {q, p}], (4) 

where in particular the constraints 

Dgg>0, Dpp>0, DggDpp-D^p-(^'f >Q ( 5 ) 


hold, which provide the necessary and sufficient conditions for the dynamics described by 
0 to be in Lindblad form and therefore completely positive pT|[38) . An alternative way 
to get the same positivity condition is to ask the generalized Heisenberg uncertainty relation 
{q^)t{p^)t — (({p, 9 })t/ 2 )^ > fi.^/4 to hold for any time and any initial state | |4^ . 

The first two assumptions are implicitly or explicitly taken in all the Markovian 
approaches. A further natural requirement is that the interaction with the environment does 
not depend on the position of the oscillator. 


Assumption 3 (Translational invariance). The reduced dynamics is invariant under 
translations apart from the harmonic potential term. This requirement is equivalent to the 
validity of the classical equations of motion for the mean values of position and momentum. 


By applying the generic translation qv^q + x, pv^pto the generator 0 we see that 
all the terms are invariant with the exclusion of the term containing the harmonic potential 
— [g, {q, p}] and the last term — 4^o-7m) qqjej-efore, the above assumption 

is satisfied if and only if kq = 7m- The same conclusion is reached by considering the 
equations of motion for position and momentum and asking them to be equivalent to the 
classical equations in which the momentum is proportional to the derivative of the position. 

The result of the first three assumptions is therefore that the Liouville operator has the 
structure 0 with kq = 7 m > 0 and > 7m/4; moreover, the constraints 0 hold. In 
particular, the Hamiltonian part of C turns out to be 

-ffm = -ffo + ^ {<7,P}, ^ 1^ 

where, besides a contribution in the form of the free Hamiltonian of a harmonic oscillator with 
a strictly positive frequency Hm, one has an additional term in the form of an anticommutator 
proportional to the damping constant. 

The evolution equations for the mean values and second moments of position and 
momentum then read; 


d(9)i _ {P)t 


d(p)t 


dt 


m 


dt 


'm^m{q)t -7m(p)t, 


(7) 











CONTENTS 


6 


d(g^)t ^ (fag})t 
dt m 

d(p^)t 


+ hD, 


qq^ 


dt 

d({g,p})t 

di 


= -"i^m({P>9})t - 2'Y^(jP)t + hDp 


( 8 ) 




- ‘2mCl'^{q^)t - 7m({g,rf)t - 


QP' 


The dynamical matrix giving the evolution of the mean values (|7]i has eigenvalues — 7111/2 and 
—7m/2±i-ym^ ~ 7m/4> which naturally leads to introduce the damped frequency Wm of the 
mechanical oscillator in terms of its bare frequency 


Wii 


= t - 




(9) 


Let us note that according to > 'y^/d > 0 we have ruled out the case = 0, which 
corresponds to a quantum Brownian particle, that is a massive particle not bounded by a 

for a recent review). 


potential in a translation invariant environment 118 41 -451 (see 

At this stage we further have to determine the diffusion coefficients Dgp, Dgg and Dpp 
appearing in Q. We will rely on the study of features of the equilibrium state, but to avoid 
the known incompatibilities with translation invariance im we formulate the equipartition 
condition in a weaker form, not asking the existence of an equilibrium Gibbs state generated 
by Ho. 

Assumption 4 (Equipartition). At equilibrium the mean kinetic energy and the mean potential 
energy have to be equal. 

Since the eigenvalues of the dynamical matrix associated to (|8J have a positive real part, 
existence of a unique attractive equilibrium state is granted, and thanks to the linearity of the 
equations the equilibrium state is actually Gaussian and determined by the mean values at 
equilibrium. Then, our equipartition condition is 


2 m 


= 2 "^^m(g^)eq, 


( 10 ) 


which gives equal weight to the mean kinetic and potential energy at equilibrium. By setting 
in the time derivatives equal to zero we come to 


Dgp — 


m7„ 


D, 


qq. 


Moreover, the second moments at equilibrium turn out to be 


(g )eq — „ ( Dgq + 


D 


pp 


(P^) 


27 . 

- ^ (n 

eq = [Dpp 


+ wfCl'^D, 


({g:P})eq = -hmDg 


( 11 ) 


( 12 ) 


99 / 


We exploit finally the residual freedom we have in the choice of the diffusion coefficients 
to get a Gibbs state as equilibrium state. However, as we already noticed, it cannot be state 
generated by Hq, and we replace it by a generic effective Hamiltonian. 

Assumption 5 (Gibbs state and temperature dependence). The equilibrium state has a Gibbs 
form which is determined by an effective Hamiltonian independent from the temperature. 

As we shall prove below, this assumption implies that the diffusion coefficients have the 
expressions 


Dqg = 


7m (2W + 1) 

2muJm ’ 


D„ = l^(2N + l). 

ZUJrn 


D nn — 


7n 


4Wn 


(2iV + l) (13) 
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and that the equilibrium state has the form 
= Tr{e-/3^^m} ■ 

In these expressions, is the mechanical Hamiltonian (j^ and N > 0 represents the mean 
number of excitations in the equilibrium state, namely 

N=^ -(15) 


Let us prove the above statements. Thanks to the Gaussianity of the equilibrium state. 
Assumption 1^ means peq oc exp|—ciT| for a suitable quadratic Hamiltonian (without 
the liner terms, because the means of position and momentum have to be zero), say H = 
^ + i rhOf + ^ {q,p}, with := — ^ > 0 in order peq to be a trace-class operator. 

Then, the eigenvalues of H have the form hui {n + 1/2) and, without changing peq, we can 
redehne c and H in such a way that Cj = Wni and c = f3, a positive constant which can be 
interpreted as the inverse temperature of the equilibrium state of the mechanical oscillator. 
Then, the mean number of excitations has the expression (0. By equating the mean values 
determined by peq with the expressions (12i, after some algebraic manipulations we get 




mrfm 

ifi^ 


mDn 


I t^pp 


2mD, 


qg 


{2N+lf = - 

7, 




^ ^ 2 
^qq^pp ^ ^qq 


The right hand side of the last equation is greater or equal to 1 by Q and (111. To have 


H independent from the temperature implies that the coefficients Dqq and are both 
proportional to 2N -f 1 with temperature independent coefficients. The equations above, 
together with Cj = Wm, give the expressions (13 1 and H = iTm. 


2.2. Master equation for the mechanical oscillator 

From the previous results we see that a central role is played by the mechanical Hamiltonian 
Hpn, which appears in the commutator part of the Liouville operator and determines the 


equilibrium state (14i. It will be very useful to diagonalize explicitly iFm by introducing 
a suitable mode operator. By dehning 

1 . „ . . UJ^ i 7m 


Om = 


(mH^g + irp), 


\/2mhuj, 

; get that the mechanical Hamiltonian ([^ can be written 


r = 


Wn 


(16) 




(17) 


The dimensionless quantity r has modulus equal to one and [ani,a)„] = 1 is satished. The 
inverse formulae giving g, p in terms of am, are 


q = 


2muj„ 


r Om + ra. 


t 


I mMTf 


p = i\ 


2UJrr 


(18) 


Let us note that the form ( [l4| of the equilibrium state does not come from a direct 
requirement, but rather it follows from all the considered assumptions. In particular we stress 
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the fact that the operator is not the Hamiltonian of the isolated oscillator, but includes a 
term containing 7 ni which comes from the interaction with the bath. Combining 0 and ( [T^ 
we have in particular 

1), ^({<z,p})eq = -|^(27V + l), (19) 


(P^)eq _ 1 ^2 / 2\ _ (nAT 


indeed satisfying 0. We see that the term related to damping gives a negative contribution 
to the equilibrium mean value (i?m)eq arising from energy exchange with the bath. The 
Lindblad operators Rj appearing in now read Ri = \/7m {N + 1) Om, R 2 = VlmN 
so that the Liouville operator can be hnally written as 


C[p] — ^ [Tfni; p] Tni(-^ “f 1) [ ^mP^i- 


- 5 ! 


1,P} 


■ ImN 




41 


it 


^p} 


( 20 ) 


Let us stress that, despite the fact that the expression p0|) has the form of the generator for 


an optical oscillator 126), the relations (|16|, (|18|l connecting Um, with g, p account for 


the description of a mechanical oscillator. A master equation for a mechanical oscillator 
with the Liouville operator (j2^ and the relation (16 1 between mode and position/momentum 
operators was already proposed in fT^ Sects. 6, 7], inside a scheme of canonical quantization 
of dissipative classical systems. The introduction of ( [T6| in order to obtain a quantum 
description of the mechanical oscillator complying with all the natural physical requirements 
is a key result of this section, which we will later exploit to consistently treat optomechanical 
systems. 


3. Langevin equations for the mechanical oscillator 

So far we have obtained a quantum master equation in Lindblad form for the mechanical 
oscillator, only relying on general physical constraints and symmetry requirements. However, 
optomechanical systems are typically dealt with making use of quantum Langevin equations, 
which provide a suitable and powerful approach for linear systems p7| ; in such a framework 
not only the system of interest appears, but also some quantum noises representing the 
environment. It is a general result that for any master equation in Lindblad form it is possible 
to introduce in a rigourous way a unitary dynamics involving the system of interest and 
suitable quantum Bose fields, which at the level of the reduced dynamics of the system exactly 
reproduces the master equation. That is, these quantum Bose helds effectively describe 
the thermal environment affecting the mechanical oscillator and the system/held dynamics 
is given by a unitary time evolution operator satisfying a quantum stochastic differential 
equation of the type introduced by Hudson and Parthasarathy 0. Within this formalism 
the Heisenberg equations for the system operators provide the quantum Langevin equations, 
while, as shown in p9) , the Heisenberg equations for the Bose helds give the input-output 
relation of Gardiner and Collet | [2^[28l . We thus obtain in a unihed framework all relevant 
physical information 0. Finally, we shall show in Section [3^ that this approach allows to 
treat also non Markovian effects and to introduce noises with non hat spectrum. 

Let us start introducing the Hudson-Parthasarathy equation or quantum stochastic 


Schrodinger equation |31|, which gives the evolution equation for the unitary dynamics 
involving the system of interest and a quantum Bose held. The proper mathematical 
formulation of this equation relies on the formalism of quantum stochastic calculus 
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For the Liouville operator p0| ) the associated Hudson-Parthasarathy equation reads (see 
e.g. or 1^ Sections 11.2.2, 11.2.7]) 


dU{t) = “ a|nd.Bth(i)) 

-^((2iV + l)al„a„, + iV)di|(7(t), (21) 

with (7(0) = 1, i7„i given by ( [T7| , and i?th(i) a Bose thermal held satisfying the canonical 
commutation relations 


[Bth{t), = min{f, s}, Bth(s)] = 0, (22) 

and the quantum ltd table 

dBth(f) dBl^{t) = {N + 1) dt, dBtV(^) dBth(t) = iVdi, 

dBth(f) dBth(f) = 0, dBth(()df = dBl^{t)dt = 0, 

with N the positive quantity introduced in ( fTS] !. The commutation rules are better understood 
by introducing the formal held densities: di?th(i) = &th(() df. Then, these densities satisfy 
the standard canonical commutation relations 


[^th {t ), (s)] =S{t- s), [6th {t), 6th(s)] = 0. (24) 


Equation is a quantum stochastic differential equation in Ito sense and the second 
line of (21 1 corresponds to the ltd correction. The solution (/(f) of (211 is a family of unitary 
operators on the overall Hilbert space which represent the dynamics of the closed system 
corresponding to “mechanical oscillator plus held”. An heuristic, but more familiar, picture 
can be obtained by using the held densities introduced above. The formal expression of the 
unitary evolution is indeed 


(/(/) = ¥ exp| J {o■mbl^^{s) - a]jj6th(s)) dsj, (25) 


where ^ denotes the time ordered product. From this formal expression one sees that U (t) 
is the time evolution operator for system and held in the interaction picture with respect to 
the free held dynamics. The thermal held Bth therefore represents the environment, say the 
phonon held interacting with the mechanical oscillator. 

It is possible to show that the physical thermal held Bth{t) does not admit a Fock 
representation. However, it is useful for computations to have a hand a mathematical 
representation of i?th(/) in terms of two commuting Bose helds Ai and A 2 in the Fock 
representation | [29]|48| . This means that such helds satisfy the canonical commutations rules 
[Ai(f), A|(s)] = ^ijmin{f, s}, [Ai{t), Aj{s)] — 0 and that there exists a common Fock 
vacuum, i.e. a normalized vector e(0) annihilated by all these operators: Ak{t)e{0) = 0 for 
k = 1,2. The held dehned by 


Bth(f) = ^/^NTiAl{t) - y^Alit), 


(26) 


satishes the canonical commutation relations ( |22| l and the Ito table ( |2^ . It is known in 
quantum held theory that there exist non unitarily equivalent representations of the canonical 
commutation relations; indeed, Bth{t) cannot be obtained by unitary transformations of Fock 
helds. 

Fet us now consider as state of the held the A-held vacuum e(0). In such a case taking 
the partial trace over the Fock space of the helds, which corresponds to take the trace over the 
environmental degrees of freedom in open quantum system theory, the reduced system state 
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is given by p{t) = Trenv {U{t) p(0) 0 |e(0))(e(0)|17(<)'^}, with p(0) the initial state of the 
oscillator. Thanks to the reduced dynamics of the mechanical oscillator can be shown to 
obey exactly the master equation ( [20| p0| . Further, we have the important relations 

(e(0)lBth(t)Blj^(s)e(0)) = (A^ + 1) min{f, s}, 

(e(0)|B^^jj(f)Bth(s)e(0)) = Nmin{t,s}, (27) 

(e(0)|Bth(f)Bth(s)e(0)) = 0. 

It is worth noticing that the thermal parameter N does not appear in the commutation rules 


of the field Btt, but rather in the quantum correlations (27 1 . This expresses the fact that N 
depends on the “state” of the field or, more precisely, N determines a non-Fock representation 
of the canonical commutation relations. Note furthermore that the vacuum e(0) is not 


annihilated by the fields Bth(t), but it plays the role of a thermal state |48 Sect. 6]; no 
vacuum state exists for a non-Fock Bose field. 


5.7. Quantum Langevin equations and input-output relations 

Relying on the previously introduced formalism we are now in the position to obtain the 
so-called quantum Langevin equations, providing the stochastic evolution for the system 
observables in the Heisenberg picture. For a generic system operator X we denote as usual 
the Heisenberg picture as X{t) = U{t)^XU{t), with U{t) the unitary operator describing 
the closed dynamics of system and environment. Differentiating this expression by the rules 


of quantum stochastic calculus, essentially summarized by the ltd table (231, one obtains the 
quantum Langevin equations for the relevant system operators, namely for the mode operator 


da^{t) = - (^iwni -b a^{t)dt - ^/^dBth{t). 
By (18 1 we get also the equivalent equations for position and momentum 

dq{t) = df -b dCq{t), 
m 

dp{t) = - {mn^q{t) -b lYnP{t)) dt -b dC'p(f), 

in which we have introduced the Hermitian quantum noises 

flln 


C,(t) = - 




(28) 

(29) 

(30) 

(31) 


Cp{t) = (i7th(f) - Blit)) , 

where r is the pase factor defined in ( [T6| ). By ( [22] l the new noises obey the commutation rules 

[Cq{t),Cp{s)] = i;i 7 mmin{f, s}, [C'g(f), Cq(s)] = [Cpit),Cp{s)] = 0. (32) 

A fundamental advantage of the considered formalism is that, thanks to the unitarity of U{t), 
the transformation X i—)■ U{t)^XU{t) preserves all the commutation rules among system 
observables, in particular the Heisenberg relations between position and momentum, as can 
be checked also directly relying on ( |32| l. Warranting preservation of these fundamental 
commutation relations is indeed a basic step in providing a true quantum description of a 


dissipative dynamics |26 Chapts. 1, 3]. 


We now consider the Heisenberg picture for the thermal fields and we define 

75^ (t) = U{t)^B,i,{t)U{t). 


( 33 ) 
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While Bth{t) represents the field before the interaction with the oscillator, the so-called 
input field, is the field after the interaction, the so-called output field. We stress in 

particular that an important consequence of the Hudson-Parthasarathy equation is the identity 
B°^^(t) = U(T)^Bthit)U{T), Vr > t, which warrants the fact that the output fields obey 
the same commutation relations as the input fields, namely ( |22) l; in other words, both the 
input and the output fields behave as free fields. By differentiating the three contributions in 
U{ty Bth{t)U{t) according to the Ito rules, one gets the input-output relation 

dB°y\t) = dBth{t) d- am(l) df. (34) 

The linearity of the Heisenberg equations of motion allows for an explicit solution 



-’th 


it) = - 


- IWn 


T + iWn 


Bth{t) 


7n 


7ir 




2 I ■ Jo 


'dBth(s) 


(36) 


The explicit expressions for q(t) and p{t) can be easily obtained from (181 and (35 1 . 


3.2. Field state and non-Markovian dynamics 


In the Markovian approximation considered so far, the temperature enters the theory only 


through the parameter N defined in (151. This approximation can be described stating that 
the system actually sees a flat noise spectrum, or more precisely the system is only affected 
by the value of the bath spectrum at the frequency Wm- A more general and physically more 
realistic situation is to allow for a structured noise spectrum and this can be achieved without 
any modification of the unitary dynamics ( [ST) and of the related Langevin equations and 
input-output relations. To this aim it is enough to change the state of the field by taking 
mixtures of coherent states |30 Let us note that considering such a mixture of coherent 
states for the description of the state of the field is actually analogous to consider a state with 
a regular P-representation in the case of discrete modes (see e.g. p^), as explained below. 
This modification is new in the context of quantum stochastic calculus and will imply that the 
reduced dynamics of the oscillator is no more Markovian, in the sense that a closed master 
equation in Lindblad form for the statistical operator cannot be obtained. 


3.2.1. The field state. In order to consider a more general field state let us first introduce the 
Weyl operators ||^[^ for the Fock A-fields, defined as 


r ^ ' ^+oo 'j 

>Va( 5 ) =exp|^y^ fffe(s)dAj'(s)-h.c.j 


with gk square integrable functions. The operator WA{g) is unitary and the property 
Afc(f)WA(5)e(0) = /g ds(7fe(s))A’yi((7)e(0) holds, so that the action of a Weyl operator on 
the Fock vacuum gives a coherent state. Therefore WA{g) is nothing but a displacement 
operator for the Bose fields | |49) . Relying on (261, we can introduce a Weyl operator also for 
the P-field 


>Vt(/) = exp 



f{s)dBl{s) 



(37) 
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where / is a locally square integrable function and T denotes a suitable large time, which 
we will let tend to infinity in the final formulae describing the quantities of direct physical 
interest. Now, yV’T(/)e(0) is not a coherent state for the Sth-field, but its relevant moments 
can be computed by using the ^-field representation ( [26] i: 

(WT(/)e(0)|i3th(f)WT(/)e(0)) = f f{r) dr, 


(WT(/)e(0)|Bth(f)Sth(s)WT(/)e(0)) = / f{u)du /(r) dr, 


(Wr(/)e(0)|St\(s)Bth(f) WT(/)e(0)) = iVmin{f,s} + / f{u)du fir) dr, 


(38) 


(Wr(/)e(0)|Sth(f)S^\(s) WT(/)e(0)) = {N + 1) min{f, s}+ /(u) du / /(r) dr. 

Jo Jo 

A crucial step is now to consider / to be a random process and to take the state of the 
field characterizing the environment to be 

= E [WT(/)|e(0))(e(0)|WT(/)^] ■ (39) 

Again, in the final formulae we will take the limit T —+oo. This is nothing but an 
analogue of the regular P-representation for the case of discrete modes. Indeed, in the case 
of a single mode the Glauber-Sudarshan P-representation of a state p |26 Section 4.4.3] 
is defined by p = f d^a P(a, a*)|a)(a|. If the pseudo-density P is allowed to become 
negative and singular, then any state can be represented in this form. When P is a true 
probability density, one speaks of a regular P-representation and the Glauber-Sudarshan 
formula describes mixtures of coherent states, including in particular thermal states |26 p. 
113]. In a probabilistic language, which is more suitable for generalizations to stochastic 
processes and fields, the fact that a state p has a regular P-representation can be rephrased by 
saying that it can be written as the expectation value p = E[|a) (a|], with a a complex random 
variable. In order to construct a thermal state it is enough to consider the case in which the 
distribution of a is Gaussian with vanishing mean E[a] = 0 and second moments E[a^] = 0, 
E[|a|^] = (7^. Then, p = E[|a)(a|] turns out to be a thermal state |26 Section 4.4.5]. 

By analogy, to construct a thermal field state with a general thermal spectrum we take 
/ to be a Gaussian stationary stochastic process with vanishing mean, E[/(f)] = 0, and 
correlation functions 

E[/W/(s)]=0, E[7(f)/(s)]=:G(t-s). (40) 

Thanks to stationarity, the function G(f) is positive definite, so that according to Bochner’s 
theorem its Fourier transform 


Giu) = 




e-''^‘G(f) dt 


(41) 


is a positive function, which we assume to be absolutely integrable, thus implying a finite 
power spectral density for the process. Since the field state denv, defined by ( |39| ), turns out 
to be Gaussian, we can characterize it through the means and the correlations of the thermal 
field Pth. which are immediately obtained from p8] ) and the properties of the process /. The 
only non zero contributions are given by 

(Pt^ij(s)Pth(i))env = ^niin{f, s} -f /" du f drG{r-u), 

Jo Jo 


{Bthit)Bl^is))enw = [N+l)min{t,s} + f du f drGir-u). 

Jo Jo 


( 42 ) 
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To better grasp the physical content of the new state and of the formulae ( |42| ) let us 
introduce a set of held modes as in | |4^ . Using a complete orthonormal set {/i„} in 
we can expand the held in terms of discrete independent temporal modes by dehning them as 

/■ + 0O 

= / 


hn{t) dBthit). 


We then obtain {ch„)em = 0, (c^^)env = 0, together with 


lim (4 Chjenv = N + E |(/ln|/)z ,2 


^+oo 


T—¥--\-oo 

where we have dehned the positive quantity 
N{iy) = N + G{iy) 

and hn{v) is the Fourier transform of hn{t)', by normalization 


hn{v) N{v)Av, 


(43) 


\hn{v)\'^dv = 1. So, 

the reduced state of the single mode Ch^ is exactly a thermal state expressed in the P- 
representation. If we take hi and /12 having non overlapping Fourier transforms we also 
get lim 7 ’_).+oo(c)ji''^ 2 )env = 0, which means that these two modes are independent. Then, 
N{i') is naturally interpreted as the mean number of phonons in a given field mode Ch well 
peaked around the value v of the frequency and held modes with different frequencies are 
independent. A value of N(v) different from zero in a neighbourhood of v implies that 
the mechanical oscillator can absorb from the bath phonons with energy around hv. On the 
contrary, the approximations are such that the oscillator can emit phonons of any frequency, 
even when Niv) = 0. The physically relevant quantity is now the combination of the two 
non negative contributions N and G(y), rather than the values of the individual quantities. 
Note that the Markovian reduced dynamics of Sectionj^can be obtained either by considering 
the non-Fock representation for the thermal field, thus assuming a strictly positive W > 0 in 
(26 1 and taking G'(;/) = 0, or equivalently by considering a standard Fock representation and 


formally taking the limit of constant spectrum G{y) in all the physical quantities. 


5.2.2. The equilibrium state of the mechanical oscillator. According to the dehnition of 
reduced dynamics, the time evolved state of the mechanical oscillator is still obtained 
by taking the partial trace with respect to the held degrees of freedom p(t) = 
limT-).+oo Trenv {U{f) (p(0) ® CTenv) However, at variance with the case in which 

the state of the held was taken to be the A-held vacuum, by taking the time derivative of this 
expression no closed evolution equation is obtained unless N{v) is constant. Not to have 
a closed time-homogeneous equation for the reduced dynamics is indeed a signature of the 
non-Markov features of such a dynamics. 

In spite of the difficulty of not having a closed master equation, the study of the reduced 
equilibrium state, namely peq = linit_>.+oo can still be afforded and its expression 
enlightens the physical role of the various parameters. Indeed, thanks to the requirement 


E[/(f) /(s)] = 0, one has that equipartition in the sense of (lOl still holds. Starting from 


the explicit forms of position and momentum in the Heisenberg picture (see ( [35] l, ( fTS] !) one 
can check that the equilibrium mean values of position and momentum remain equal to zero, 
while the variances are still of the form ( [T9] l with N replaced by the effective mean number 
of excitations 

N[v) 


NeS = 


7m 

27r 


■dv. 


Notice that if the quantity N{v) introduced in (43 1 is 
corresponding to the Markovian case, then A^eff = N. 


(44) 

taken to be the constant N, 
This result suggests that the 
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final Markov approximation should be valid when G{v) is approximately constant in a 
neighbourhood of uj-^. In fact equation (44 1 represents a smearing of N(y) around the 


frequency of the mechanical oscillator Wm, the more peaked the smaller the damping 
constant 7„i- Non-Markovian effects can only be relevant if G{v) appreciably varies in a 
neighbourhood of width 7 „i around Wm, being suppressed with decreasing 7i„- 

Since the equilibrium state is necessarily Gaussian, by comparing ( |44) i with ( [T4) i we get 
that the new equilibrium state is again a Gibbs state with respect to the same Hamiltonian H^, 
but with an effective inverse temperature /3eff defined by setting — 1)“^. 


3.3. Properties of the quantum noises and quantum stochastic Newton equation 

Let us now come back to the quantum Langevin equations for the position and momentum 
operators, so as to better understand their physical meaning and the role of the noises. In 


order to study the properties of the noises we transform the Langevin equations (29 1 , (301 in 
the form of a stochastic Newton equation. 


To this aim we first have to consider the quantum noises (31 1 appearing in these quantum 
Langevin equations. The commutation relations ( |3^ for these noises, which are state 
independent, guarantee the preservation of the canonical Heisenberg commutation relations. 
Their quantum correlations do instead reflect the physical properties of the field state cTenv and 
can be obtained starting from the H-correlations ( |42] i. Note that Langevin equations for a 
mechanical oscillator of the same form and with two noises obeying the same commutations 
rules d3^ were used also in however, the two point correlations used in p^ , taken 
from |17| , arise from approximations in the Caldeira-Legget model which are valid only for 
medium/high temperatures, while in the present treatment they are deduced from the state of 
the phonon environment and are valid at any temperature. 

We stress the fact that in the present formulation the momentum operator is not related 
to the time derivative of the position operator according to the classical relation, but rather 


through (291 where the quantum noise C'g(f) explicitly appears. However, the connection to 
the classical formulation is not completely lost. In fact from ( |2^ we can derive the relation 

gfe) - q{tl) _ 1 /■*" p(f) ^ Cgfe) - Gq{tl) 

t 2 - ft t 2 - h m t 2 -tl 

By ( pT) and ( |42) i, the mean value of the r.h.s. of the equation above vanishes, while its variance 
is given by 


{{G,{t2) - G,{h))% 


^7n 


mWni (t2 —ti)\2 


1 


r+oo 2 ^sin 

-oo {*2 - h) 


N{iy) du 


(t2 — ti) 

SO that in particular also the variance goes to zero for growing t 2 — G. Then the quantity 
v{t) = p{t)/m can actually be interpreted as the “coarse grained” velocity of the mechanical 
oscillator. 


If we use the formal field densities bth{t), with commutation rules (24i, take as 


starting point the quantum Langevin equations ( |29) l and ( p0) l and eliminate the momentum, 
we can rewrite the quantum Langevin equations in the Newton form: 

mq{t) + + mn^q{t) =C{t), (45) 

where we have introduced the formally Hermitian quantum noise ^(f) 

^(t) = Cp{t) + m"fmCq{t) + mCq{f). 


(46) 
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Most importantly the commutation relations for this noise take the singular expression 

= 2imnj^^6{t - s). (47) 

While the expectation value of this noise with respect to the field state is zero, its 
symmetrized correlation function can be computed from the relations 




dtds 


(C,(f)C,(s))e„v = 


1 




(Cp(f)Cp(s))e 


dtds 

TV + i — s) + Re G{t — s) 


(48) 


1 

^7m 


d^ 

dtds 


({C,(<),Cp(s)})e„v = - 


UJn 


+ - ) <5(f — s) + Re G{t — s) 


— ImG'(t — s) 
7m 


(49) 


and has the expression 

^(UW,CW})env=^ 

Z Win 

+ 2m/i7i, 




d^ 


dtds^} 


N + - ] 6(t — s) + Ke G{t — s) 


-ImG(f-.) 


(50) 


Note that (45 i and ( |47] i were already introduced in |26 Sect. 3.1.2] and |37|, where the 
commutation rules ( |47| i were actually enforced by the requirement of preservation of the 
commutation rules between position and momentum. However, at variance with previous 
approaches, here we have provided an explicit construction of the quantum noise ^{t) in terms 
of a quantum Bose field, based on a rigorous mathematical construction. 


We stress the fact that the stochastic Newton equation (451 is mathematically purely 


formal due to the presence in (46 1 of Cq{t), which contains the formal derivative 6th(f) and its 


adjoint. Moreover, if one were to take ( [45] l, ( [47| and (50 1 as starting point for the construction 
of the quantum Langevin equations for position and momentum, then one should complete 
(45 1 , which is an equation for q{t) only, with a suitable definition of p{t). The standard 


choice in this respect, considered for instance in |[T 25 is to take p(t) = mq[t). This 
works out fine as far as the commutation relations of position and momentum are concerned. 


However, in this case the equation of motion (45 1 and the structure of the noise ^{t) obeying 


(46 1 imply that q(t) contains singular quantum fluctuations, so that it is not a well defined 


operator. Also p is then not a well defined operator and its variance is actually infinite. Then, 
one has to regularize the momentum, by subtracting the noise responsible of this divergence; 
this is what our construction does. The identification of the momentum is given implicitly 


through the first canonical equation (29 1 , which corresponds to the coarse grained velocity, as 


discussed above. No divergency appears because the whole construction is based on the well 
defined unitary evolution (|2T]l. 


3.3.1. Consistency of the quantum noises. It is important to stress that if a set of quantum 
Langevin equations is considered as starting point for the description of a stochastic quantum 
dynamics, commutations rules and symmetrized correlations of the noises cannot be given 
arbitrarily. In particular, independently of the considered system, if {‘Ci(f)} is a set of 
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operator valued noises, the quantum correlation function (fO)env has to be positive 

definite fST), in the sense that 


E 


^+oo 


^+oo 


df / dit' 


( 51 ) 


for every choice of the “smooth” test functions {hi}. Since we can always write 

^ + I 0 (^ 0 ] , ( 52 ) 

the necessary positivity condition introduced above becomes a consistency condition between 
commutation rules and symmetrized correlations. 


Relying on (48 1 , (491, as well as the commutation relations 


for the noises Cq and 


Cp, one can immediately check this fact for the model at hand. Also for the singular noise ^ 


constrained by (461 one can show that the expression is positive dehnite. These 


results are due to the fact that the noise helds have here been explicitly constructed in terms 
of the quantum Bose helds, so that commutation rules and correlations are not postulated, but 
rather follow from the mathematical expression of the model. 


3.3.2. The noise correlations. 
correlation of the noise ^ by 

= I I 


For the model at hand we denote the Fourier transform of the 


^+oo 


dfe (U(f+ s),^(s)})e 


I that according to ( [4^ and ( [50| it reads 

mh'y, 


R{u) = 


2uJr, 


^ + (wm + 


N{v) + - 


+ —)■ —v), 


(53) 


(54) 


where {v —> —v) means to add the same contribution with v replaced by —v. Note in 
particular that R{v) is an even function of the frequency. 

In our treatment, which gives rise to the expression (54i, the interaction with the 


environment is described in terms of exchange of quanta with the bosonic held representing 
the phonons, see ( |25] l. In models in which the system of interest is coupled to other 
harmonic oscillators, by some approximations it is possible to arrive to a quantum stochastic 
Newton equation like ( |45| ), but with a different noise spectrum. A reference expression often 
considered in the literature 126 (3.3.9)], 137 for the quantity R is given by 


Rgz{v) = TThJ{iy) coth 


J(^) := — k{v)v. 


(55) 


The quantity k{v) contains information on both coupling constant and density of modes of 
the bath in a neighbourhood of the frequency v, J{v) is often called spectral density. Also 
in the case of this choice of the noise correlations, it is possible to show that the equilibrium 
mean of qltY' diverges and therefore the identihcation of the momentum with mq{t) is not 
possible, but some regularization is needed. A typical choice in this context is k{v) = 7 „i |26 
(3.1.1)], this is equivalent to J{v) oc v, which is known as Ohmic spectral 

density. The function k{i') must be even due to the dehnition (531 and stationarity. Also the 
commutations rules ( [47| and the positivity requirement (51 1 still have to hold, leading to the 
requirement Rgz{v) — mhy^ |z/| > 0, satished at any temperature by taking k{v) > 7 in > 0 


for all V. This requirement tells us that in order to have a consistent model satisfying (45 i, at 
least at large times, and preserving Heisenberg commutation relations one cannot consider a 
spectral density with gaps inside the expression (55 1 of Rgz- 
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By a suitable choice of N{v), it is possible to get R{v) very similar to Rgz{v), apart 
from the low temperature limit. For instance, by taking 

(56) 


N{v) = \v\ k{v) 


+ 1/2) 7ni - 1) 


we obtain 


Riv) = Rgz{v) + mh 


7m 

2u}„ 




(57) 


note that the difference is independent from the temperature. In ( |57) l the compatibility with 
the commutation relations is guaranteed by the first term in the square brackets; so, in ( |5^ 
there is no restriction on the choice of k{v), apart from k{v) > 0, and even spectral gaps can 
be introduced. 

Besides the case (56 1 , the freedom in the choice of N{v) in (54i allows to model 
quite different environments, for instance with a sub-Ohmic or super-Ohmic spectral density 
| |22|j24t , or with a structured occupation spectrum. In optomechanical systems the quantity 
R{v) enters experirnental quantities, as in the cases of homodyne/heterodyne detection 
discussed in Section 


4.3 


so, in principle it is possible to test the form of R{iy). However, 
essentially a frequency window of width 7 ni around Wm is experimentally relevant. Some 
results at room temperature ||52) seem to indicate a non-Ohmic spectral density around a 


very interesting possibility, but not enough to discriminate between RgzUA ^nd the form (571 
for Riv). At zero temperature, corresponding to N{v) = 0 in our case (54i and to /3 —>■ +C !0 
in the difference is most evident, namely 


R{v) = mh'j^ 


^1 + ^^ 

2Wm 


versus 


Rgz{v) = mhk{v) \v\ 


> 7n 


(58) 


(59) 


These expressions are very different, but to discriminate between them one has to consider 
very low temperatures and a ratio 7m/wm that is not too small. 


4. Cooling and emission spectra of an optomechanical system 

As an application of the quantum description of a mechanical oscillator developed so far 
we consider the simplest optomechanical system |[^|^|^|6 25 namely the mechanical 


oscillator is a mirror mounted on a cantilever and coupled to the light in an optical cavity by 
radiation pressure. The cavity is of high quality, without thermal dissipation other than the 
one due to the coupling between cantilever and phonons and tuned in such a way that only 
one electromagnetic mode is relevant. Strong laser light is injected and some light is allowed 
to leave the cavity so that its spectrum can be analysed. 


4.1. The optomechanical model 

The micro-mechanical oscillator (the mirror) is described by the operators q, p as in ( fTS] ! and 
by the Hamiltonian (j^. The cavity mode is described by the operators Oc, aj and by the 
free Hamiltonian liwcajac- The free electromagnetic field is in a coherent state describing a 
perfectly monochromatic laser of frequency Wq; however we use the equivalent description of 
inserting a source term for the cavity mode in the Hamiltonian and of taking the external field 
in the vacuum. The final optomechanical Hamiltonian takes the form 

Honi{t) = TTm + hwcaluc — hgoqalac + ihE — ace'^”*) . ( 60 ) 
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Note the trilinear term giving the interaction between the position of the mirror and the number 
operator of the photons in the cavity, which represents the radiation pressure interaction; the 
coupling constant is usually expressed as go = uJcjL, where L is the length of the cavity. The 
laser power is P = ftuoE'^/"fc, where 7 c is the cavity decay rate and E the laser amplitude. 

In order to include the cavity mode interacting through radiation pressure with the 
mechanical oscillator, as well as the emission and absorption of the light from the free 
electromagnetic field, the Hudson-Parthasarathy equation ( [2l] i is modified as follows; 

dU{t) = I acdBj^{t) - h.c.^ 


7m 

2 


((27V + 1) + N) dt- Y alocdt >U{t). 


( 61 ) 


Here Bth is the thermal field satisfying ( [2^ , while B^m is an independent Bose field in the 
Fock representation, describing the electromagnetic field outside the cavity. The relevant Ito 
rule is dBe^{t)dBj^{t) = dt, while all the other possible products vanish. Now U{t) is the 
unitary operator describing the dynamics of the two interacting oscillators and the fields. The 
latter are in a factorized state given by the tensor product of the thermal environment state 
([39)l and the electromagnetic vacuum; 


denv — ^env C) |eem(0)) (Cem(O) | ■ 


(62) 


It is convenient to eliminate the laser frequency working in the rotating frame and 
introducing the unitary operator V{t) = which upon differentiation obeys an 

equation of the form ( [hi) albeit with iTom(f) substituted by 

iTin + hAoalac — Hgoqalac + ihE (aj. — Uc) , (63) 

with Aq = uJc — ojQ the nominal detuning. For a generic system operator X we define 
X{t) = V{tyxv{t), so that by differentiating according to the rules of quantum stochastic 
calculus, as done in Section [XT] we get the following quantum Langevin equations 

dac(f) = (^iAo + y) ac{t) + igoq{t)ac{t) - iE^dt - (64) 


(65) 


dCp(f) 


of Section 3.1 


we 


as well as 

dq{t) = df + dCq{t), 
m 

Mt) = {-mfl^qit) - 7mP(t) + hgoal{t)ac{t)) dt 

where Cq and Cp are given by ( jTTj l. Defining the output fields as in 
have besides (34 1 the input-output relation for the electromagnetic field 

dB°^(t) = dBem(f) + ^/^e~'‘^°^ac(t)dt. (66) 

In the case of a very intense laser, that is E^ large, the dynamics can be linearized in 
a neighbourhood of the equilibrium mean values, determined by self-consistency from the 
means of the linearized form of the quantum Langevin equations. The equilibrium mean 
value of the momentum is zero, while setting = {ac{t))eq, we find 


c = - 


iE 

¥+Ta’ 


(„) = ^30 ICI' 


(67) 


where we have introduced the effective detuning A, 


A = Ao - go(g)eq = Wc - go{q)eq “ Wq. 


(68) 
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By inserting the equations ( |67| ) into ( |68l l we obtain the self-consistency condition 
{A-Ao)(^+ aA + hg^E^ = 0; 


(69) 


this cubic equation determines A as a function of the laser parameters Ag and E. 

In writing and solving the linearized quantum Langevin equations it is useful to have 
dimensionless and selfadjoint system operators. It is therefore convenient to set 


Xit) = 


CaJ(f) -f Cac(i) 


- (g)eq), p{t) 

-72 Id, Y{t) = 


pit) 


i(C4W -Cac(f)) 


(70) 


(71) 


72|d ' ” 72|C| 

Then, the linearized quantum Langevin equations turn out to be 

dw{t) = Aw{t)dt - dQ{t), w{t) = (g(f), p{t), X{t), Y{t))'^ , (72) 

where the superscript T means transposition and the dynamical matrix is given by 

A A 
A A 

-Aimc -iT-r 


A = 


A — 

■^mr. — 


Am — 


0 rim 


Ar.= 


0 


(73) 

G^/ Wm/flm 

oj’ 

-A - 


(74) 


The quantity G, having the dimension of a frequency, will play the role of effective coupling 
constant and is given by 


G = 9o ICI 


2h 

mCCni 


(75) 


(77) 


so that in particular it depends on the effective detuning A through given in (67 1 . The vector 
of noises is given by the following field combinations: 

Qiit) = +h.c., Q,{t) = -ph.c., (76) 

Q,{t) = +h.c., 

Q 4 (t) = +h.c., 

where t is the phase factor defined in (16 1 and the quadratures Qiit) and ( 52 (f), apart from a 
multiplicative factor due to the change of dimensions, coincide with the noises 

Note the different structure of the two dynamical sub-matrices in ( [74| . Indeed the 
former describes a mechanical oscillator and the latter an optical mode, corresponding to 
different interactions as discussed in Section The same choice is taken, for instance, 
in |[T]|2l|3|2^|^[^, but not in ||5]|3[T^. 

The linearization around the equilibrium state is meaningful provided one can ensure 
the existence of such a state. Its stability conditions can be obtained by applying the Routh- 
Hurwitz criterion to the equations for the mean values, which correspond to the system (721 
with the noise term dQ{t) suppressed; the result is the couple of conditions 

,2 


G^UJmA < If' 


75 

4 


-A' 


(78) 
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for A > 0, and 




|A| < 


7c 7n 


7c + 7m 


+ + A=) + 


(n.i - ^ - A*)' 

7m +7c 


(79) 


for A < 0; there is no restriction for A = 0. The same stability conditions have been found 
in 125), as their equations for the mean values agree with ours. 


4.2. Energy fluctuations and laser cooling 

To introduce the fluctuation spectra of position and momentum of the mechanical oscillator 
we use a formulation tailored for (classical or quantum) processes starting at time zero and 
we define the gated Fourier transforms 0 




Ai/t 


dBflt), 


Vt Jo 

for the Bose fields as well as for the relevant system operators 


i = th, em, 


(80) 


= 7 t h ' 


ivt 


Wi(t)dt, i = 1, 2, 3,4. 


(81) 


Here T is a large time going to infinity in the final formulae to recover a stationary situation. 
Then, the spectra of the fluctuations of position and momentum of the mechanical oscillator 
are defined, in analogy with the classical case |[54), by the quantum expectations 


1 


Sqiv) = Fi(T;-u)}), 

1 —^ 

Spiv) = lim h{F2{T;u), F2(T;-i/)}), 


T —^-j-cxD 2 


1 


S„iu)= ^\un^-{{F,iT-u), F2iT--v)} + {FflT--u), F^iT'o)}). 


(82) 


(83) 


Let us stress that, while useful, these definitions do not correspond to some continuous 
monitoring of position and momentum, even though Sq{v) is directly related to the observed 
optical spectra as we shall see in Section [43| 

The Fourier transformed equations of motion corresponding to ( |72) can be solved by 
purely algebraic manipulations and the vector F{T-,iy) can be computed; due to the length of 
the expressions the result is reported in Appendix A| To compute the spectra above we need 
also the held correlations, which we give in ( |A.9| |. 

Due to the vanishing of the held cross-correlations, the spectra (82 1 , ( [8^ decompose in 
a thermal and a radiation pressure contribution according to 


Sqiu) = Sl^iu) + Sq,iu) = 

2 


Spiv) = ^Sl^iu) + Sl^iu). 


(84) 


By inserting the expressions ( |A.l| i, ( |A.2| i into the definitions (82 1 , (831 and by using (A.9 1 we 
get, by some computations, the expressions for the spectra of the fluctuations: 


S7iv) = 


Sira^mG 

‘2\div)\^ 


A'^ + A + A 

4 


(85) 
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2Wm^^m \d{v)Y 


-G^ 


hm \d{u)\ 

I (^m +v(^ 


■ Tn 

V\LOr„-\ — 


N{i 


2)+(- 


\ + {y + ^f 




( 86 ) 


(87) 


=-^ (^)+ 




2 \d{u)Y 


+ ^ - ) - ly^ciu^m + l^) 


+ (i/ -)■ -v) 


( 88 ) 


where (12 —> —:/) means to add the same contribution with v replaced by — ^ and the quantity 
R{v) is the Fourier transform of the quantum correlations of the noise given in (54i. The 
denominator d{v) is the characteristic polynomial of the dynamical matrix A given in (73i 
and ( |74| l: 

d(i/) = det (^ + i^l) = + i - A + i - G^WmA. (89) 


Note that the quantities (85 1 -( 90 1 are non-negative as they should be in order to give a 
sensible decomposition of the spectra. Another useful way to write 8 ^( 1 /) is by putting in 
evidence its difference from S'‘*'(z^); the resulting expression is 


- Sl\u) = (Ni,.) + 1 


rim \d{l^)\ 


1 G^A + 

2 2UJrri 


j.- _ - Al 

4 


+ iiy^-iy)y 


(90) 


4.2.1. The peaks in the fluctuation spectra. A relevant role in determining the properties 
of the system is played by the denominator d{v) (89i; indeed, the quantity ^^[A^ — 
(z/+ i7c/2)^]/d(z/) is sometimes interpreted as the effective mechanical susceptibility |25 
Eq. (17)]. Most importantly note that the zeros of d{v) determine the positions and the widths 
of the peaks of the fluctuation spectra; even though in principle they can be obtained by 
solving the fourth order algebraic equation d{v) = 0, it is much more convenient to have 


simple expressions, even if approximate. An analysis of these zeros is given in Appendix A.l 
in the case in which d(y) can be written in the form 


d{v) = 


+ i “ Aggf 



I • Tm \ m 2 

+ 1 ~ ^eff 


(91) 


The stability conditions (78 1 , guarantee the strict positivity of the effective damping 
constants Tc and Tm. The quantity ui^ is known as optical spring rigidity, while the ratio 
(Tni — 7m)/7m is Called co-operativity ||^|^. 

An exact expression for the zeros is found when A = Wm, which allows us to put into 
evidence a crossing of the frequencies of the hybridized optical and mechanical modes Q. If 
also the condition 


4G2 < ( 7 , - 7^)2 


(92) 
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holds, the result is 


r. = 


7c +7n 


+ e\ 2v? - 2a; 2 + 


(7c - 7m)^ 


r., = 


7c +7n 


- e\ 2v?- - 2 lu^ + 


(7c - 7m)^ 


A 2 _ , 2 _ 




(7c - 7m)^ 

32 


where 


u^ = 




(7c - 7m)^ 
16 


- G^a;^, 


e = 


+1 

-1 


if 7c > 7n 
if 7c < 7n 


Always for A = a;m, under the conditions 




< G^u:i < uj^ + 


(7c - 7m)' 


16 


w,; > 


we get instead the result 


r = r = 

-L r -L m 


7c +7n 


AefI = 


(7c - 7m)^ 

16 




_ 2 (7c - 7m)^ , 

x± i Ul^ 


G2 - 


(7c - 7m)^ 


(93) 

(94) 

(95) 

(96) 

(97) 

(98) 

(99) 


The two alternatives in (98 1 are completely equivalent; there is no reason to associate the 
frequency to the cavity and to the mechanical oscillator, or viceversa. A striking 
feature of the case A = ui^ is the change in the behaviour of the zeros at the critical point 
7c = 7cj solution of = ( 7 c — 7m)^/4; recall that G^ depends on 7 c due to (67 1 and (75 i. 


In the general case, an approximate expression can be obtained under the conditions 

lx(A)|<l, |x(A)| 


^«1, 

7c 


7c 




4a; ^ 


< 1 , ( 100 ) 


where 


X(A) = 


G^uj^N 


+ (A - o;^)^) (+ (A + 


( 101 ) 


The result for the damping constants is 

Tin 7m + X(A) (7c - 7m) ^ 


Tc 7c - x(A) (7c - 7m) • (102) 


The expressions for and A Jj are then obtained by inserting Tm and Tc in the equations 
( A.12| i. The compatibility conditions ( A.13| l have to hold. As one can see, when A > 0 
{red detuning), we have an increasing of the mechanical damping constant, Tm > 7m, and a 
decreasing of the spring rigidity, o;^ < 
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4.2.2. The mean values at equilibrium. By integrating in their frequency dependence 
the fluctuation spectra one obtains the second moments of position and momentum in the 
equilibrium state: 


(9^)eq-« = 


2TTmCL„ 


Sq{v)dv, 


(/)eq = 


mMl„ 

2 tt 




(103) 


1 

2 


({9>ri)eq = 


h 


Sqp{v)dv. 


(104) 


All these quantities are finite since the integrands behave as j/ ^ for \i/\ +oo. Moreover, 
the re duce d equ ilibrium state of the mechanical oscillator is a Gaussian state characterized by 


(103 I and (104 1 together with (g)eq = %o ICI (p)eq = 0. 


On the contrary the integral of v^Sq{v), which would give the fluctuations at equilibrium 
of q, does not exist. This fact is related to the features of the noise in the thermal 

part and, as already noticed right before Section [3.3.1| this noticeably implies that the standard 
identification of mq with momentum is not possible. The expression of Sq{i') coincides with 
the one given in |T 25 , where however Rqz{v) with Ohmic spectral density appears instead 

still q^ does not have a finite mean and 


of R{v). While in the case of 01 g Sq{u) 


also in this case the identification of momentum and velocity is not possible. Notice that the 
expressions for Sp{v) and Sqp{v) have not been obtained before. In particular Sqp{v) ^ 0 
implies that the fluctuations of position and momentum are actually correlated. 

The mean energy of the harmonic oscillator at equilibrium takes the form 


{Hra)eq — ^ 'tnftj^{q) 


(105) 


(106) 


where the contribution due to fluctuations is given by 

= ^ [ du d- Sp{iy)) +7m'S'gp(l/)] . 

Jr 

It is convenient and natural to split this contribution into three distinct terms, distinguishing a 
radiation pressure term from the rest and further dividing the thermal contributions into two, 
putting into evidence a contribution which is not proportional to position fluctuations and does 
not have a definite sign. We thus introduce the dimensionless quantities 


as well as 


Mth(A) = 




2cJrn^n 






UJ- 


27r 


dz/ 


G27„,A 


2M(^)r 


G'^N-v 


7c7n 


,2 \ “I 




N{y) 


' im 


(uJm + v) 

+ -v) 


(107) 


(108) 


so that the fluctuation contribution can be written as 


(iT)fl —+ A/th + 7Vlth(A)); (109) 

by construction we have A/’rp+A/th+Adth(A) > 1/2. As it appears, the mean energy density 
cannot be obtained from the knowledge of Sq alone, but extra terms are present. Moreover, the 
contribution proportional to Adth(A) can be negative. Depending on the parameter values, 
the extra terms can be actually quite small. It is important to stress that the given expression 
for the mean energy of the resonator holds for any temperature of the phonon bath, including 
the case of zero temperature. 
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We further stress that there is not strict energy equipartition. This can be expected since 
the mechanical oscillator is coupled to the cavity through its position and also the counter¬ 
rotating terms contribute to the final result. In the thermal part the lack of equipartition is due 
to the terms proportional to A, which are present in S^{i/) and not in In the radiation 

pressure part the term with 17^ comes from the position and the one with comes from the 
momentum and give different contributions to the mean energy. 


4.2.3. Vanishing effective detuning. For a vanishing effective detuning A = 0 all the 
computations can be performed analytically. The second thermal contribution A4th(A) 
vanishes and the coupling constant takes the value = Shg^E'^/ For the spectra 

of the fluctuations the explicit expressions reduce to 


S7i^) = 


Tc 


(z/ — Wm) iy + Ulyff) + 




^m7n 


N{u) + i 


+ (z/ -)► -v) 


( 110 ) 


leading upon integration to 


i mnl ^ (2Aeff + 1) + 


hWynG'^ (27m + 7c) 


87n 


(7n 


4 


+ w. 


0 ’ 


1 

2m 


(F^)eq = 


hnl 


(2Weff + 1) + 


hWnyG^Jc 


87n 


(7m-|-7c) 


■ W' 


^({g,rf)eq = -^(2iVeff + l), 

4 SWm 

with Agff as in (|4^. These expressions show that equipartition of the mean energy is not 
valid just due to the radiation pressure contributions. However equipartition approximately 
holds for 7 c ^ 27 m, which is the case typically considered in many theoretical studies and 
experiments. We further have for the fluctuation contributions to the mean energy 


Afyp — 


(7m + 7c) 



(7m+7c)^ 

4 



A/jh = A^eff + 2> Adth(A) = 0. 


The mean equilibrium energy of the mechanical oscillator is thus increased due to the 
interaction with the cavity as a consequence of the presence of the strong laser in resonance. 
For the values considered in Figurej^we have A/(p ~ 1.6 x 10^ corresponding to a temperature 
of about 7.9 K. 


4.2.4. Laser cooling. As discussed in many papers |[T|[3j[5l[6l[5^, an important effect which 
can be described by the quantum models of cavity optomechanics is the laser cooling of the 
mechanical resonator. Since, as already discussed, we cannot expect equipartition of the mean 
mechanical energy, we cannot speak of temperature in a strict sense. A natural way to speak 
about laser cooling is the comparison of the mean energy of the fluctuations of the mechanical 
oscillator in the presence or the absence of the stimulating laser (corresponding to = 0). So, 
we have to study the value of the fluctuation contribution ( |109| l and to compare it to its value 
for C = 0, which is given by = (^^m)eq|^^o = ^^m {N^s + \). 

To obtain explicit analytical formulae for the mean energy we consider the case of a 
constant noise spectrum, that is N(y) = const = To actually perform the calculations 
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we need the expressions of the zeros of d{v)', here we consider the generic case given by the 
Ansatz By lengthy computations the integrals over i/ can be exactly performed, leading 
to involved formulae explicitly given in Appendix A.2 In order to describe cooling effects 
the relevant contributions can be written in the form 


Mh = ^ Q 

-L m 


Nee + - 


Mth(A) = ^/C 

-L m 


iVeff+2 


( 111 ) 


C = 


(Q + /C). 


( 112 ) 


where the quantities Q and K, are given in equations ( |A.16| l and ( A.17| i. The expression for 
Kp is given in (|A. 15). Note that, while Q is always positive, depending on the values of the 
parameters the quantity K, can be also negative. For a large choice of the parameters Q turns 
out to be close to 1 . 

In the following figures we describe the effective cooling of the mechanical oscillator, 
by considering as a hgure of merit the cooling factor 

In 
Tn. 

We study two cases, corresponding to the parameter regions for which an exact or approximate 
analytic evaluation of the different contributions to the mean energy has been provided. In 
both cases mass and bare frequency of the mechanic oscillator are taken to be m = 2.5 x 10“^° 
kg and = 27r x 10^ Hz, while the mechanical damping factor is 7 m = 27r x 10^ Hz. We 
consider a cavity of length 5 x 10“"^ m and resonance frequency Wc = 27rc/(1064 x 10“®) 
Hz, driven by a laser with a power of 5 x 10“® W. For the sake of comparison the values of 
the hxed parameters are taken from p5| . 

We start by studying the case A = ujm, studied in Section 4.2.1 in which the location 


of the poles can be evaluated exactly, provided one distinguishes two regions according to the 
value of the ratio ( 7 c— 7 m)^/ 4 G®. No approximation is taken in the expression of the integrals 
giving the mean energy. If this ratio is above one, verihed for a cavity damping 7 (. > 7 c, 
where the critical damping % is introduced in the comments after ( [99| and corresponds for 
the considered parameters to 7 c — 4.1 x 10^, the effective damping rates Fm (93 1 and Fc (94) 
are actually distinct, while Agff = and their expression is given by (95 1 . By numerical 
computations we see that the cooling factor C is a monotonic increasing function of the cavity 
damping rate 7 c, and around the starting point % the cooling factor takes the value 2.9 x 10“®. 
In the complementary region, corresponding to ( 7 c — 7 m)^/ 4 G® below one, the cooling factor 
is a decreasing function of the cavity damping rate, so that the optimal cooling is obtained for 
7 c = 7 c. In this region, corresponding to 7 c < %, we have Fm = Fc with value given in 


_ ]l, while the effective frequencies Aes and ujfg are given by the expressions (98 1 . To assess 

the relevance of the various contributions in (| 1 12 |i we report the values for 7 c = 7 c: we have 
7 m/rm 3.05xl0“^ Q ~ .997 and/C ~ -4.18 x 10“^; then, (flT^ gives C - 2.91x10“®, 
which is a very strong cooling factor. 

Instead, in Figure[T]we consider the case 7 c ^ uJm, that is a cavity damping much bigger 
than the mechanical oscillator frequency. In the exact formulae for the integrals we use the 
approximate expressions for Fm and Fc given in ( |102| l, relying on the conditions ( |100| l. The 
stationary value of the energy of the mechanical system has a marked dependence on the 
effective detuning A and the optimal cooling region, corresponding to C of the order of 10“®, 
is obtained for A < 7 c. In this parameter region Af^p can be neglected with respect to CNgg, 
unless the phonon bath is below 1 K, so that indeed the quantity C given in ( |1 12| i properly 
describes the cooling effect. When the detuning A goes to zero the cooling factor rapidly 
increases in agreement with the discussion in Section 4.2.3 showing the presence of heating 
at A = 0; in this parameter region, the cooling effect disappears also when A grows. 























CONTENTS 


26 



Figure 1. Plot of the cooling factor C for the case in which the cavity damping is much bigger 
than the mechanical oscillator frequency. We explore the dependence of the cooling factor on 
both the effective detuning A and the cavity damping rate 7c, both expressed in Hz. It appears 
that the best cooling factor is of the order 10“^ and con'esponds to ^ 5; 7c- 


4.3. Optical spectra 

We consider now the monitoring of the emitted light by balanced homodyne and heterodyne 
detection | |5^ Sect. 7.2]. The aim is to see which kind of information on the mechanical 
oscillator can be obtained by detection of the emitted light. 


4.3.1. Homodyne spectrum. The case of a perfect coherent monochromatic local oscillator 
of frequency wq with detection of the whole emitted light 149 551 corresponds to the 
continuous measurement of a field quadrature of the type 

rt 


=ie-‘V"’'s« f +h.c.-, 


(113) 


Jo 

i? is a free parameter which depends on the optical path and determines the observed 
quadrature. As a consequence of the definition we have that r?), (5(s; t?)] = 0. 
By the properties of the output fields, discussed after Eq. ( [3^ , the commutation rules 
are preserved; this gives that the output current := U{t)^Q{t]d)U{t) satisfies 

t?), t?)] = 0. This is the key property expressing the fact that 

can be measured with continuity in time. Similarly to ( [SOl l we introduce the gated Fourier 
transforms 


= ^ r?). (114) 

From the above relations we obtain the second key relation which guarantees the presence of 
commuting observables and therefore the consistency of the theory: 




( 115 ) 
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The homodyne spectrum is then given by the expression 
S{v,-d) = lim Ti 

T—H-oo 




(116) 


where the environmental state is given by ( |62] l and po is any initial state for the mechanical 
oscillator and the cavity mode. Note that this expression is nothing but the spectrum of the 
classical stochastic process representing the output, and not an ad-hoc quantum definition | 


Sect. 4]. The commutation property (115 1 implies that the homodyne spectrum S{u; -d) is an 
even function of u. 

As shown in Appendix B.l the homodyne spectrum has both an elastic and an inelastic 
component 


S{l'] d) = d), 

which turn out to have the expressions 

Sei{iy;d) = 87r7c |C|^ {s'mdf S{u), 

5inel(j^; I?) = d) + Srp{u-, d), 


with 





cosd + A 

>2 2 
sinz?) +(z^cosz?) 


[ll + iA-uf) 

(^ + (A + i.)') 



(117) 


(118) 


(119) 


Srpiv;^) = 1 + 


27cWmG^ cosr^ + Asin^?)^-I-(^cosz?)^ 




+ (A — + (A + 

' — A^^ sin2z? — A ( 7 c cos 2d — 2iv) 

d{v) + (A - + (A + 


+ Ic^mG Re 


X I ^ + A^ - i7cz/ 


( 120 ) 


Note that all the contributions are indeed positive as shown in [Appendix B.l| It is important 
to stress that the connection between Sq{i') and S'inei(^^; is far from simple. In particular 
the last contribution in ( 120 i comes from the interference of the electromagnetic part of the 
signal with the shot noise, as detailed in [Appendix B.l This is a completely new term, in 
principle detectable in experiments at very low temperatures. 

Let us further stress that different quadratures are incompatible and actually one can 
prove the general inequality p^[55) 

S'inel(j^; 'd)Sinelil^; d ± 7r/2) > 1, (121) 


which is just a form of the Heisenberg-Robertson uncertainty relations coming from the 
canonical commutation relations of the involved Bose fields. As a result quite different 
physical information can be extracted from the different quadratures. 
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The case A = 0. The first striking example of strong dependence on is in the case A = 0. 
Fori? = 7r/2 we get 7r/2) = 327ri?^(5(i^)/7c and 7r/2) = 1: only the shot noise 

contributes to the inelastic spectrum. 

On the contrary, for the quadrature with i9 = 0 we get 7r/2) = 0 and 


0 ) — 1 + ' 






Sg{u), 


( 122 ) 


where Sq{u) is now explicitly given by (1101. An important point is that in this case the 


interference term vanishes exactly and we have a direct connection of the homodyne spectrum 
with the fluctuation spectrum of the position of the mirror. This result has been found also 


but with the substitution R{h') Rgz{v) in the expression (861 for . . As 


at least in principle, when A = 0 the homodyne observation of the quadrature with 
i? = 7r/2 can give direct experimental information on the correct expression for R{v). At zero 


temperature one could experimentally discriminate between our result (581 and the standard 
proposal ( [59] l. 

In other cases the interference term does not vanish, but it can be negligible at high 
temperatures. For instance, when the interference term is negligible, at least in the region 
where N{y) ^ 1, we recover for the result given in |[^ Sect. 3]. At high 

temperatures the inelastic homodyne spectrum allows to reconstruct the fluctuation spectrum 
of position, while no direct information on the fluctuation of the momentum and on the cross¬ 
correlation is obtained. Moreover, at high temperatures we have also Sq{v) ~ by 


using the explicit expressions of (861 and R{i') (54i we get 

( 


‘5'inel(^j0) — TcTin^ 


7c 

4 


^ + (o 


4 ) 


JV(i-) + - 


+ {l' ^ -v). 


This expression highlights the dependence of the homodyne spectrum on the thermal spectrum 


N[i>) and the characteristic polynomial d{v) (891 of the dynamical matrix (731 of the full 
optomechanical system. 


Squeezing. An important information about the non classical nature of the light generated by 
optomechanical systems can be obtained considering the quadrature with i? = —7r/4. In the 
simple case of vanishing detuning A = 0 and vanishing temperature N{v) = 0, it is possible 
to show from (117 1 -( 120 1 that we have 5'inei(0; —7r/4) < 1, at least in a certain region of the 


parameters. This means that in a neighbourhood of i/ = 0 we have S'inei(*^; ~7’‘/4) < 1 and 
the emitted light is squeezed. This result shows that such a kind of optomechanical systems 
can generate non classical light az)- Note that, if light squeezing is present for certain values 
of the parameters, then the inequality ( |121| i implies that the complementary quadrature is anti¬ 
squeezed. Of course, experimentally it could be difficult to tune the values of the various free 
parameters in order to have squeezing; moreover, the elastic peak in the spectrum tends to 
hide the squeezing around = 0 in the inelastic spectrum. 


4.3.2. Heterodyne spectrum. In the case of heterodyne detection the local oscillator and the 
stimulating light are produced by different laser sources; now, the stimulating laser frequency 
ujQ and the local oscillator frequency, say p, are in general different. Moreover, the phase 
difference cannot be maintained stable and this erases some interference terms. It can be 
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shown Sect. 3.5] that the balanced heterodyne detection scheme corresponds to the 

measurement in continuous time of the observables 

I{n]t)= [ V^e-"(‘-")/2eV«+i“dSem(s)+h.c., (123) 

where a is a phase depending on the optical paths and x: > 0, represents the 

detector response function. As we shall see, the heterodyne spectrum does not depend on a. 
In the Heisenberg description the observables become the “output current” 

= +h.c. 

By the definition of I{n;t) and the properties of U{t) we get [/out(p-; ^out(/i; s)] = 0, 
which says that the output current at time t and the current at time s are compatible 
observables. 

While in the homodyne scheme the spectrum of the output is analysed, in the heterodyne 
scheme it is usual to register only the output power as a function of the frequency /i of the local 
oscillator. The mean output power of the detection apparatus at large times is proportional to 

P(/i) = lim ^ [ df Tr {/out(Ai;i)^Po ® deny} ; (124) 

the limit is in the sense of the distributions in /r. As a function of p, P(/i) is known as power 
spectrum. Note that to change p. means to change the frequency of the local oscillator, that 
is to change the measuring apparatus. In general /out(M; 0 ^^d /out(At^ s) do not commute, 
even for t = s. Then, there is no reason for the power spectrum to have some symmetry in p. 
The heterodyne power spectrum can be decomposed in an elastic and an inelastic part 

P{p) = Sel(Ai) + Sinel(Ai), (125) 


Eei(Ai)= lim / df [2Re fce‘“ / e-^ 

T->-+oo T Jo L V do / 


H 2 


- ^2 —^47r7c|C| d{p-uJo), 


^ + {t — ^o) 


where 


(126) 


1 r'^ 

Sinel(M) = lim — df Tr{/i„el(At;f)^Po ® deny} , (127) 

T^+oo T Jq ^ 


+ y|^e'l'^-“«l*+‘^(y(s)-iA:(s))ds) +h.c. 

The inelastic part of the spectrum is computed in [Appendix B.2 Again it is possible to 
identify a radiation pressure contribution and a thermal part 


Sinel(A^) — ^rp(A^) T ^th(P')- 
For simplicity we give only the expressions for k ^O: 


Srp(/r) — 1 + 


SIP {p - Wo) 


- Im 


7cWmG^ If -i{p-uJo + A) 




(^7l + (^^-ujo- A)") d{p - wo) ^ + i (m - wo - A) 


(128) 

(129) 
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- Wq) 

~ 1^0 ~ ^) ^ 


( 130 ) 


Both contributions are positive as it follows from the expressions ( 130| l and (B.7 1 . Note the 
presence of the interference term in p29[). By simple computations one can check that 


Sinel(j^ + Wo) + Sinei(wo - v) = 5'i„ei(r^; t?) + 5'inel(j^; ^ + 7r/2); (131) 


this is a fundamental relation | [5^ Eq. (9.61)] connecting heterodyne and homodyne spectra. 
Moreover, by inserting the definitions of the relevant quantities given in ( |5^ , ( |85| ) and ( [86] l, 
an explicit expression for Einei can be obtained from which it is apparent that Einei(Ar) > 1: 
in the heterodyne detection the phase dependencies are lost and it is impossible to detect 
squeezing in the emitted light. 



lo'^ 


Figure 2. Plot of the inelastic heterodyne spectmm Sjnei as a function of v for a range of 
values of the cavity damping 7 c around the critical value 7 c discussed in Section [4. 2.4| It 
appears how the two distinct peaks of the spectrum coalesce at critical value. The spectrum is 
plotted for A = a;m, while the other parameters are as in Section[4.2.4| 


As in the homodyne case, the interference term in (1291 is negligible when ^ 1 and 
we get 

7cWi„G'^ 


^inel(A^) — 1 “t“ 


Clr 


+ (/r - Wo - 


Sqiy^l - Wo). 


(132) 


When this approximation holds, the inelastic heterodyne spectrum too allows to reconstruct 
the asymptotic dynamics of the mirror through the position fluctuations. 

To explore the behaviour of the spectrum we take N(y) as given by (561 with a Ohmic 
spectral density. Then, by using the explicit expressions of S'JP and and by setting 
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v = fj, — luq, we get 
Sinel(^^ + Wo) = 1 + 


TcG" 


2\d{i^)\ 
2 


7c<G2 


■7n 


f+ (^ + A)^ 


r ^.2 


7m I . , ,, ^2 4a;„,|i.| 


(133) 


From this expression we see that the main features of the spectrum wi ll be d etermined by the 

for A = Wm we 


4.2.1 


zeros of the denominator \d{v)\ ; for instance, as discussed in Section 
can have one or two resonance frequencies depending on the value of the cavity decay rate 
7 c. In Figurej^we show this phenomenon: the two distinct peaks coalesce as 7 c increases. 
For these values of the parameters one can check that the main contribution to the inelastic 
heterodyne spectrum comes from the thermal part Eth- It can be checked that in this parameter 
region the behaviour of the inelastic homodyne spectrum S'inei(j^;0) is very close to the 
heterodyne one as depicted in Figure Let us notice that the behaviour shown in Figure 
does not uncover the whole rich structure of the spectrum which appears by exploring other 
parameter regions. 


5. Summary and outlook 

In this article we have shown how to give a fully quantum description of a dissipative 
mechanical oscillator. The combined use of master equations and quantum Langevin 
equations allows for the construction of a dissipative dynamics respecting symmetries and 
physical constraints, such as the energy equipartition at equilibrium, and subject to dissipation 
with an arbitrary noise spectrum. A crucial feature allowing for these results is that for a 
mechanical oscillator the definition of the creation and annihilation operators and in 
terms of position and momentum is not the usual one, but rather depends on the damping 
constant as discussed in Sectionthe standard result is only recovered for a vanishing 
damping constant as can be seen from (|l6|. Moreover, the quantum Langevin equations for 
the system, and the input-output relations for the noises, for both the mechanical oscillator 
and for the optomechanical system, given in Section 3.1 and Section [4T| respectively, need 
not be postulated: they are nothing but the Heisenberg equations of motion determined by 
the Hudson-Parthasarathy unitary evolutions pT) and ( [6T] ). In this framework it appears that, 
in order to preserve the Heisenberg uncertainty relations, the momentum operator can be 
interpreted as the time derivative of the position operator only in a “coarse grained” picture. 
An help in comparing our approach to others and in discussing the structure of the noises 


comes from the quantum Langevin equation in Newton form, see Section 3.3 which at the 
price of introducing singular noises does not contain the momentum operator. Indeed in the 
quantum case important constraints on the correlation functions of the operator noises come 
from the fact that they need to be positive definite and compatible with the commutation rules 
of such noises. In this formalism, we are further able to introduce a field analog of the P- 
representation for the state of the environment and this opens the possibility of treating an 


arbitrary noise spectrum as done in Section 3.2 


Our description of the mechanical oscillator is not very different from other proposals 
at medium and high temperatures of the phonon bath. Differences become relevant for 
very small temperatures. Indeed the dynamics we have constructed is fully “quantum” at 
all temperatures and this opens the possibility of constructing models of optomechanical 
systems which are reliable also in a deep quantum regime. As an example we have studied 
a prototypical system: a mechanical resonator interacting via radiation pressure with a single 
optical mode in a cavity. For this case we have given explicit general formulae for the 
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fluctuation spectra of position and momentum of the mechanical resonator and for the mean 
mechanical energy at equilibrium. By using detection theory in continuous time, we have 
obtained the full expressions of the homodyne and heterodyne spectra of the emitted light. For 
not too low temperatures, usual results are recovered, such as laser cooling and connection 
between the light spectra and the fluctuations of position of the mechanical component. 
However, our description is valid also at very low temperatures, when semi-classical reasoning 
is not valid and the observation of the spectra of the emitted light is not giving a direct 
measurement of the mechanical fluctuations. 

Many generalizations are possible |57 -|^|, which could benefit of a systematic and 
consistent treatment. The simplest generalization is to include imperfections in the detection 
scheme and noise in the stimulating laser light ||^|^[^[^. But also direct detection can 
be included pO) or the entanglement between resonator and optical mode can be studied. 
Moreover, the whole theory has in some sense “modular” properties and can be applied to 
more complicated systems, say when more mechanical resonators and more optical modes 
are involved. 
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Appendix A. 


Computation of the fluctuations 


The Fourier transformed equations of motion corresponding to (72i can be solved by purely 
algebraic manipulations; essentially the problem reduces to compute the inverse of the matrix 


A + iv\. By using the characteristic polynomial d{v) of the dynamical matrix A (891 the final 
result for large T turns out to be 

, 2 


Fi(T; j.)- 
F2{T- j .)- 




d{v) 


20Jn 




/u;m7c 


d{v) V 2 
ZjP(r;z.), 


F4(T; v ) 





Zl^{T-u), (A.l) 

(A.2) 
(A.3) 


(A.4) 

(A.5) 


ZT(T-v) = 




+ (I - ^ ‘ t)) - 


7n 


LOm — 1 I I — G WmA 


(A.6) 
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Z7{T\v) = (i (i/ - A) - y) - v)^ 


+ (i + A) - y ^ 1^ + Wo), 

nT i.. ,.^t , /^ A , ; 7c\ C 


(A.7) 


Zl^{T- ^) = (a - - i I) ^ 4^^(wo - ^)t + (a + ^. + i I) ^ 44(^ + wo). (A.8) 

To compute the spectra ( [8^ , ( [8^ we need also the field correlations. From the 
correlations (|4^ and the fact that _Bem is a Fock field in the vacuum state we get 

(44(^.)44(i.)t)e„v = 1, {Bl{v)^Bl{v)U = N{y), 

(44(444(i/))e„v = 0, {Bl{u)Bl,{u)^U = N{u) + l, 

while the cross-correlations involving both and Sem vanish. Also the fluctuations of the 
cavity mode operators and the correlations oscillator/mode could be computed by using all 


the components (A. 1 1 -( A.41, but we do not study them in this work. 


Appendix A.l. The zeros of diy) 

Let us assume that d{v) has two zeros of the form = uj^ — irm/2 and = Aos — irc/2 
with Wgg f 0 and Aes f 0; by the property d{v) = d{—V), the other two zeros are —V^ and 
—Fc. Therefore, we can write d{v) in the form (91 1 or 

d{v) = [v - i^m) + Fm) [u - Uc) {v + TC) ■ (A. 10) 

By equating this expression to ([89)t we get the algebraic system 


PniH- Yq — '7 cH“ Tm; 

Tc Wmf + Fjn Wcf = + 7n 


A' 


r r — Q' 

-L m i- r — J 


4 


2 I 7c 
4 

2 

' ' 7c7m, 


(A.ll) 




A" + ^ - G^oj^A. 


The stability conditions ( [78l l, ( [7^ guarantee Fm > 0 and Fc > 0. By assuming Fc f Fm, 
from this system we get in particular 




r — r 

-L r i- n 


m 2 _ 7m 2 _ 

^eff p _ p 

-L r -L rn 


Be - Fn_ 
7c - Tc 
F. - F„ 


A^ - (7e - Fc) 


4 

Be - 7n 


(A. 12) 


An exact expression for Fm and Fc can be found when A = w^. We 


The case A = a. 

study only the case of d{i/) of the form (A. 101 with four distinct zeros. 


In the case Fc Fm we set x = Fc — Fi„ and insert (A. 12 1 and Fc + Fm = 7c + 7m into 
the last equation of the system ( |A.ll| l; in such a way we get 

-f [l6Wj^ - (7c - 7m)^] + 64G^W,^ - 16w^(7c - 7m)^ = 0. 


Then, by using the solution of the equation for x"^ and Eqs. (A.121, we find the result (931- 


(96i. By imposing Fc, Fm, A^^ to be real and strictly positive and Fc f Fm, we get the 


necessary and sufficient condition (92 1 . 

By the choice Fc = Fm, from the system ( A.l l| i we get directly the result (98 i, (99 1 , 
together with the conditions (97 1 . 
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An approximate expression. To compute approximately Tm we adapt a suggestion given 
in l[^[^ and based on an approximation of the mechanical susceptibility. In the expression 


of taken from (89 1 we make the approximation (r-m + A + i^) (i/m — A + ~ 


'2 


(u;m + A + i '^° — a + i) and we solve d(i/m) — 0 for under the 

conditions (1001, (|101|l. By using also the first equation of the system (A.ll 1 we get the 


expression (102 1 for the damping constants. Once we have Tm and Tc, we can compute uj^ 


eff 


and Ajf from the equations (A. 121, which do not contain approximations. 


For consistency we need the positivity of Ajj and which means the positivity of 
the right hand sides of equations (|A.12[). Under the approximations (1001-( 102 1 , this gives 


A2 


X(A) 


(7c - 7m) 


CJ. 




X(A) 


^ + (1-x(A))(1-2x(A)) 


(7c - 7m)" 


OJ. 




(A.13) 


<1-X(A); 


because x(A) has the same sign as A, conditions (|A.13|| give true restrictions only for 
A > 0. We see also that conditions ( |A.13| l are violated for A positive and small. In this 
situation the cavity is overdamped and the decomposition of d(iy) takes the form d(iy) = 
(ly - uj^ + i %^) (^ + + i (^ + i (i^ + i ^); we do not study this case. 


Appendix A.2. Computation of the mean mechanical energy 

The integrals over v in ( |107| l, (1081 can be performed by the residue method, under the Ansatz 
(91 1 and N{i') = Aeff. First we set 

2 \ 2 

= ... ■ , ... .... 

(A. 14) 

7c" 


L± = 


u,V = ^er + r^/4. 


With this notation we have 

G^7 c [ G^tu^A 




- A^ ± Ai 


7m^ni T 7c ( A + ^ ) + 7m7c (7m + 7c) 


+ l(7c + 7m)L 


(A.15) 


where |i/c| Wm\ is given by the last of (A.ll i. The thermal contributions A/th and Adth(A) 
are given in (111 1 in terms of the expressions 

L + 2d^ 


Q = 


“ '(7c+7m) 




+ 


2T,D^ 


Wc.f kml" 


16 

.,,2 


+ 27c^^m + 27m ( A^ + ^ 


7c ( A' + ^ + 7c7m (7c + 7m) 


(A. 16) 


K. = 


G^A r(^" + w) (T-^m) ■ 
2r,D^ [ 


‘2utm kcl" kml" 


7m^2jjj + 7c ( A + j + 7c7m (7c + 7m) 






7cIA2 + 2^ 


■7i 


m V 2 


■7c 


2 ' / 2a;ni 

Note that, while Q is always positive, JC can also take on negative values. 


(A. 17) 
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Appendix B. 

Appendix B.l. 


Computation of the optical spectra 

The homodyne spectrum 


The homodyne spectrum (116i involves the quantity (114i; by the rules of 

quantum stochastic calculus we can compute i?), which turns out to contain the 

quantities ( |A.3| l, ( |A.4| l. Then, by integration we obtain 

gr(^;^ V> ICI Sint? e-^/2^ (g 

uy 1 

Qt = EtY,{v\'d)T BJ^{v) + Eth{-v]'d)TB'^i^{-v)\ 

= -Eemy, t?) 

Eth{l^; t?) = -G-v/TmTc (^ + ^ (^ + ^m)) 


Eem{v,'d) = - ^ 


^ - i (i/ - A) iwni7cG'^e‘’’L(j/; t?) 


L(i/;t?) = 


^ + i - A) 

Asint? + + ii^) cost? 


^ + i (i/ - A) 


(B.2) 

(B.3) 

(B.4) 


Note that L{—v; t?) = L(t/; t?). The key relation (115 i together with [Bf (i/), Bj = 1 
implies 

[Qy (z/; t?) + Q^iv, t?), Qr (-1^; t?) + Q'^{-v, t?)] = 0, 
which is equivalent to 

|T;th(z^;t?)|"- |T;th(-t^;t?)|" + |T;em(i^;t?)|'- |£;em(-i^;t?)|" = 0. (B.5) 

By long computations this relation can be verified also explicitly by using the expressions of 

T;th(z^; 1?) and i?). 

By using ( |B.l| i, ( |B.5| l a nd (|A.9|l, fro m (116 1 we get the decomposition of the homodyne 
spectrum expressed by Eqs. ( |117| , ( |ll8] l with 

5th(z/;t?) = |T;th(^;t?)|' (^iV(z^) + + |T;th(-^;t?)|' [ n {- v ) + 0 , 

(B.6) 




Note that 5'th(t^;t?) > 0 and Srp{i';'d) > 0. To compute the thermal part we note that 
\Ethiv-, z?)|^ can be written by using Riu) (54i; by taking from ( 861 , we getEq. (119 1 . 

To compute the radiation pressure component of the spectrum, we need the square 
modulus of Eejn (0, which is the sum of two terms. So, we have the square modulus 
of the first term (the shot noise), the square modulus of the second term (the signal) and the 
double product (the interference term): 

iE.„Mr = i+ 


7+ ("-*) 


+ ^niTcG^ Re 


ie-2i’’ (^-i(z.-A))+i(f-i(z. + A)) 
d{L'){f+i{v-A)) 
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By inserting this expression into (|B.6|l we get 


= 1 + 


U! 






Re 


(:^ + (A-.f)(^ + (A + .f) 

ie-2ii9 _ i (jy _ A))^ + i — ii^)^ + A^^ 

2d{iy) (^ + (A-i/)^) 


Asinr? + + ii^) cost^ 


d{iy) 


+ {iy^ -v) 


Finally, by elaborating the argument of the real part and by using the expression (851 for 
S'JP(i/) we get Eq. (120 1 . 


Appendix B.2. The heterodyne spectrum 

By a procedure similar to the one used in Appendix A| and Appendix B. 1 in the limit of >f i 0, 
xt —)• +(X), we get 

T+Kf-wo-A) i%o^7c|CI^ 


/inei(i^; t) ~ e‘“V^ 


g-f (t-s)+i/is 


'y — i (p, — u!o — A) md{p — tuo) 


'2 ^0 + A) ^ i%0^7cC g-2i(/x-LJo)s-2ia 

^ — i{p — ujQ — A) md{uJo — p) 


d-Sem ( 5 ) 


le 


1 filmic 

C 

V 2mWm 

_(i(wo - p) 


— + i (/i — Wo + A 


c 


X (^ + i (M - CCO - Wn,)) e-2KM-o)-2ia 
X -i{p-uJo + A)^ - i(F - + Wm)) dBth(s)| +h.c.. 


d{p - Wo) 


By using the field correlations ( |A.9| l we can compute the heterodyne spectrum; again, by 
the vanishing of the field cross-correlations, the thermal contribution and the electromagnetic 
contributions decouple in the expression of the spectrum. By some long manipulations and 
by recalling that the limit in \\21\ is in the sense of distributions, we get Eq. ( |128| l with the 
thermal part given by (130 1 and 

iwm7cG^ T +A)|^ 


Srp(/i) — 


1 + 


2d{p — Wo) ^ -f i (/r — Wo — A) 


, ,2„, 2^4 

4|d(/r-wo)|^ 


(B.7) 


which becomes \\29\ by expanding the absolute value and using (|110|i. 
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